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$ Г. ВВЕДЕНИЕ 

Одной из актуальных задач биотехнологии является совершенст- 
вование существующих [1] и разработка новых, в том числе мемб- 
ранных [2] методов селективного извлечения целевых продуктов био- 
синтеза из культуральной жидкости. В последнее время в научной 
литературе был рассмотрен новый метод интенсификации селективно- 
го транспорта заряженных компонентов [3-й] _‚ в частности, белков 
[6-7] через синтетические полимерные мембраны, состоящий в спе- 
циальном подборе химического состава растворов и создании перепа- 
дов концентраций микроионов на мембране. Для полного выявления 
возможностей этого метода необходимо создание его математической 
модели [8] , [9], в частности, описание переноса заряженных бел- 
ковых макромолекул при совместном действии градиентов концентра- 
ции сразу многих растворенных веществ, то есть диффузиофореза в 
многокомпонентном растворе. Наиболее простое решение этой задачи 
может быть получено, если вычислить градиенты электрохимических 
потенциалов всех ионов в растворе и затем рассматривать диффузио- 
форетическое движение частицы как алгебраическую сумму движений 
под действием ионов отдельных сортов - ионоборетических движений. 
Математический аппарат для такого рассмотрения развит в [10-12]. 
Там вычислены ионофоретические подвижности частицы с тонким двой- 
ным электрическим слоем (ДЭС) при произвольной величине поверх- 
ностного потенциала. Однако для описания переноса белков более 
адекватным представляется рассмотрение другого. подмножества ха-.. 
рактерных параметров задачи, а именно случая слабозаряженных, час- 
тиц с размером, соизмеримым с толщиной ДЭС. Так, при описании 
электрофореза белков традиционно используется формула Гейри [13] 
В данной работе будет получен аналог формулы Генри для ионофоре-- 
тической и дийфузиофоретической подвижностей. 


} 2. ПОСТАНОВКА ЗАДАЧИ 


Исходные уравйения для описания ионофореза под действием 
ионов сорта © (#=1,..,и) выведены в [11,12] ; они включают урав- 
нения Стокса для ползущего течения жидкости (при малых числах 
Рейнольдса) и уравнение непрерывности для потока ионов с учётом 
их избытка в ДЭС: 


м Ге, ЧТ СУ: -%/р.)] = 0, 
Чу = 0, (т) 

2 | -2. 
1 0-0 7 - Яр+Ё, = О, ф = с; (е "Р1у9. 
Здесь П - динамическая вязкость жидкости, р. - коэффициент 
диффузии ионов сорта © , 2, - их алгебраическая валентность 
(положительная для катионов и отрицательная для анионов), Ж- 
электрохимический потенциал, Е. = объёмная концентрация, а} - 
скорость течения жидкости относительно частицы, Р - потоковое 
давление, #. - потоковая ре #. = равновесный потенциал 
ДЭС, измеренный в единицах Ртл на (он считается заданным). 


Граничные условия на поверхности непроводящей непроницаемовй для 
жидкости частицы суть: 


1-а Эх вы 


(а - радиус частицы), то есть условия прилипания для скорости и 
условие отсутствия нормального потока ионов. Граничные условия 
на бесконечности зависят от того, какое конкретное явление мы 
рассматриваем; для задачи ионофореза необходимо задать однород- 
ный градиент электрохимического потенциала ионов @-го сорта 
(индекс сорта, если не оговорено противное, здесь и далее опус- 
каем): 


9 

Е = Вс»б (3) 
91 1+ е> 

(здесь @ - угол между осью симметрии системы и радиусом - век- 

тором Т ‚ направленным в данную точку из центра частицы) и 

условие отсутствия полной силы, действующей на частицу: 


Е-0. (4) 


Поле скорости должно быть регулярно на бесконечности; однород- 
ная на бесконечности скорость движения жидкости относительно 
частицы (взятая с отрицательным знаком) и будет при этом скоро- 
стью ионофореза. 


$ 3. ПРИБЛИЖЕННЫЙ МЕТОД РЕШЕНИЯ 


Система уравнений (Т) в математическом отношении представ- 
ляет собой систему линейных дифференциальных уравнений с пере- 
менными коэффициентами; получить ее общее решение, справедливое 
при всех значениях параметров, не представляется возможным. 
Будем разыскивать приближённое решение, рассматривая в качестве 
малого параметра безразмерный потенциал ДЭС ф .- В теораи элек- 
трофореза такой подход разработан в статьях Генри [14], Овер- 
бека [15] , Буса [16] , а в теории диффузиофореза и ионофореза 
он ранее не применялся. Представим все полевые величины в (Т) 

в виде рядов по целым неотрицательным степеням этой малой вели- 
чины, не выделяя её явно: 


= - + 1, + 
Р = Р» + Р +Р2+.., (5) 


> > > — 
ф южн ++ 
А = № + | + М: +... 
Тогда для величин Е. первого и второго порядков получим 
соответственно: 
Г.) фе = = 


и 17. пл, о 

119 ей т. - Яр = = 0, (6) 
ды = -Ч\[аф (р, - ,/О)], 
су 17 = 0, (®) 


я: - УР м - 0. # = - СФУ, 
ди, = @и[-= Чр(Ярь- 55/Ру2 + 2. - %/Р)], 
@йу 5 1, = 0, (8) 
12 то! 5 3, - Ур, +Ё -0; = - сефбм са фи, 2. 
Наша цель - вычислить 17’ с точностью до членов второго поряд- 
ка; при этом надо учесть, что в силу граничных условий члены ну- 


левого порядка в скорости и давлении отсутствуют. Таким образом, 


нам необходимо вычислить поля Ве й [а р 17. й 9, . 


Поле Ко даётся выражением: 


№ = Вс (1 + 2) (9) 


(в формулах (9) - (15) сферический радиус Т измеряется в едини- 
цах радиуса частицы а). Для вычисления остальных полей необхо- 
димо решать неоднородные уравнения с постоянными коэффициентами. 
Боли ограничиться полями дипольной симметрии (которые только и 
могут дать однородную на бесконечности скорость), то нетрудно 
получить общие решения этих уравнений методом вариации постоян- 
ных. Общее решение уравнения Пуассона для электрохимического по- 
тенциала 


ди = 92) соб 
имеет вид: 


и = с 0 (СА; л + Аь) - 1 - 90, 
КА 


й 


то есть суммы общего решения однородного уравнения и частного 
решения неоднородного, причем последнее выбрано так, что экспо- 
ненциально убывает за пределами ДЭС. Применение этой формулы к 
первому уравнению (7) дает: 


р 2 
и, = <» 0Ва [+240 - 5), 99 
1 


где Ре 
Т = { фсх) +: са, # (х> = 2х2 + х* 
Общее Е неоднородных уравнений Стокса 
41$ = 0, , 
1 3оЁ оф 7 — УР = рес), фе СЫ мп.0, 0) 


имеет аналогичную структуру [17] : 


2 
у, = = 0СА, - 252+ ьт . № + И), 


о осА- № АРА 


> © 
с 
> 
>< 
+ № 
^ 


где частные решения вида 


р т -. а? ил В кА р |. 
У. се = 75 ое Сх- = ты 


м. 
У-= 4х [+9 (-х+ № 2, —)+ (2) 


40 ох 2. 
‚2 2 
* Те 69 (х-55 а, 


[@ (1) = Зах со (2 + х)+ осо ( 2 № 2=) 


Хх 12 


также экспоненциально убывают за пределами ДЭС. 

Преобразуем (ТТ), (12) так, чтобы непосредственно получить 
скорость ионофореза. Для этого учтём, что для регулярности ско- 
рости на бесконечности необходимо положить Аз = 0; в силу ус- 
ловия (4) А„ также равно нулю (поскольку Р = - ДХТаи Ау 
[17] >. коэффициенты Ат и Ао определяются из условий (2): 


А, - 24, + У, 4-0; -А,-А, +, 4)=0, 


откуда получаем, что интересующий нас коэффициент Ат, то есть 
модуль скорости на бесконечности, равен 


Мы = А, - 54-5, 4. 


Подставляя сюда выражения (12), приходим к формуле: 


У = Эт ) е]х [+ (СЗх+1 + 2х) +, (3+ -4*°)], (13) 


которая даёт вклад в скорость ионофореза того порядка по потен- 
циалу ДЭС, какого порядка член в объёмной силе подставлен в 


неё. Применяя эту формулу к уравнениям (7) и (8), получаем со- 
ответственно: 


4х фо у <), (14) 


ах фо [ 1+, (<) + 


+ ооо + . яч 442 1, (х,4)], (15) 


х 


ше Г 00 =- Ох 4/2 х* + 1/2х*, 


$ 


ро Зоб-хиех8, Ч бон ххх + 
+1746 ро = бхЧ-юиуй + Зу-к ух. 


Нетрудно показать, что с увеличением толщины двойного слоя 
Мл, убывает и при бесконечной толщине ДЭС обращается в 
ноль. 


$ 4. СТРОЕНИЕ РАВНОВЕСНОГО ДВОЙНОГО ЭЛЕКТРИЧЕСКОГО СЛОЯ 


Чтобы практически воспользоваться формулами (14), (15) для 
расчета ионофоретических подвижностей, необходимо подставить в 
них выражения для равновесного потенциала ДЭС слабозаряженной 
сферической частицы с двойным слоем произвольной толщины. Это 
выражение должно включать члены вплоть до квадрата `5 -потен- 
циала включительно (5 -потенциалом будем называть потенциал 
поверхности частицы, измеренный в единицах Ю ИЗ )›. если ли- 
нейный по $ -потенциалу член не зависит от состава электро- 
лита (зависит только от его ионной силы), то квадратичный по - 
потенциалу член зависит и от состава раствора. Чтобы вычислить 
эти члены, запишем уравнение Пуассона-Больцмана 


(т6) 


2 с2 г (17) 
+ 


ео Фрая 0. (18) 


Воспользовавшись малостью ъ -потенциала, разложим (17) в ряд 
по степеням я : 


2 2 
о (19) 
2 2 а 

где 2“=АЛ9Я с ЕЮТ- квадрат обратного дебаевского радиуса, 
а величину ^ = 3/04 са можно назвать асимметрией раствора. 
Величина А обычно по абсолютной величине не превосходит Т и 
только для бинарного 3:Т электролита достигает 2. Разложим те- 
перь у в ряд по степеням $ -потенциала: 


р = $ +4, +... 


и получим для отдельных членов разложения уравнения: 


Аф, = =“, , (20) 


Дф, = =“ (ф, - ^ф,./2). (21) 


Величина Ф, должна удовлетворять краевым условиям (18), а. 
-нулевым краевым условиям на обеих границах. Решения уравнений 
(20), (21) имеют вид: 

м жа „-С“-аа) 
фо = хе ь (22) 


-х х 


9, хх = Те -][ хе +Гюе, (23) 
1х 1 х 2 х 


где >< 
Т, =[ еа)-е *“], ста); 1. 0= 44 леорче Я, 


а переменные Х и У измерены в единицах 22”. Теперь для по- 
лучения ионофоретической подвижности не в виде ряда по степеням 
‚ ав виде ряда по степеням параметра \х нужно перегруп- 
пировать члены в разложении (14), (15). Подстановка $, в (14) 
даст единственный член, линейный по < ; подотановка Ф, в (14) 
и ЧФ. в (15) даст члены, квадратичные по | Подотавлять. 
$. в (15) не нужно, так как это даст члены Ос® у, кото- 
рые не учитывались в расчете. 


$ 5. РАСЧЕТ СКОРОСТИ ДИФФУЗИОФОРЕЗА В МНОГОКОМПОНЕНТНОМ 
РАСТВОРЕ. 


Скорость движения частицы под действием градиентов концен- 
трации одновременно всех растворенных веществ равна сумме ско- 
ростей ионофоретических движений, то есть членов вида (14) и 
(15), просуммированных по всем сортам ионов. Интегралы, входя- 
щие в эти выражения (с учетом (22), (23)) зависят от 924 и могут 
быть вычислены и протабулированы с помощью ЭВМ. Удобнее, одна- 
ко, предварительно провести указанное выше суммирование, ис- 
пользуя для Г, конкретные представления 


В. = ЮТ( ус, / Я ‚В 
В; = (Ус; /С; +2 Уф), В. = Ва, 
и лишь затем табулировать возникающие интегральные выражения. 
При этом скорость частицы будет иметь традиционный вид [18] 
суммы двух вкладов - вклада, обусловленного электрофорезом в 
поле диффузионного потенциала микроионов, и вклада за счет ос- 
мотической компоненты диффузиофореза. Первый вклад для слабо- 
заряженных частиц с ДЭС произвольной толщины уже рассчитан‘ в 
теории Овербека - Буса - Вироема [19] . Расчет второго вклада 
до сих пор был проведен только для частиц с тонким ЭС [20,21] у 
Суммирование: коэффициентов: при интегралах в_(14), (15) по 
всем сортам _цонов дает (после замены В; на В; ): 
и. ас; 2, В. ___8 ВТ Е (веа)? 
=] 31 6% й а 2 ’ 


п а*е. 228. д.2 ЮТ = (жа в [а (ва) 


2 У 9 Т 2. 
Здесь Е =- ЕТУ Ра Я = напряженность поля диффу- 


зионного потенциала микроионов, Т= <, 2; — ионная 

* 
сила раствора. Скорость движения. частиц с точностью до членов 
порядка (/( $) включительно запишется в виде 


и: Е ГЕ 8 ЮГЕ 
УЕ Зе 


й. 52 (57) > жа), (24) 


10 


где 
(кот “/х _-сх-ка) 
фсеа)=- 9” хе з +. (а), (25) 
_ а ‚ © ак ›-о-га = и -(4- а) 
# (2а)= # (о) у ке ее? х 


— сх 


-ха) х >, —(ч-2а) 
ее Чью 


= -(х-ма) 
$ (ка) -- фто +В | ке” фе" 


„о (27) 
«Есо- СЕ Зка)- Басре Ес} Еее, 


Подчеркнём, что выражение (24) справедливо в случае много- 
компонентного раствора произвольного состава. 

Первый член суммы (24) представляет собой скорость электро- 
фореза в приближении Генри [14] , второй член обусловлен ва- 
лентной асимметрией электролита и возникает в теории электрофо- 
реза Овербека [15] . Относительно функции в третьем чле- 
не суммы ранее было известно, что при 32а -> <> она должна 
обращаться вТ [20], а при достаточно больших ха , согласно 
[21], 4 -4-24/2 ха + ОГ (ау. ] 

Вычисление функций фСео), 92а) и $, С2а) может быть 
проведено следующим образом. Все интегралы, входящие в выраже- 
ния (25)-(27), заменой переменных приводятся к виду 


©оесо 


< к -© +4 

` в" Ю, соч жабо (16° ЗК ое 4, (28) 
00 

где ^;, ;=12-рациональные функции своих аргументов. Интегралы 

вида (28) вычисляются приближенно с помощью квадратурных формул 

наивысшей алгебраической степени точности 


те” Ро ах = р К, К, хо, 


феи ррооракау = 22. К, К ВЕ оь во, 


31 


где координаты Х; 


узлов интегрирования и значения коэффи- 
циентов ; приведены в [22]. 


В таблице Т представлены результаты расчетов для т = 15; 
для сравнения там же приведены значения выражения 6 о. + (жа) 
из монограйии [13] ‚ рассчитанные М.Гориным в Т94Т году. 


Таблица 1 


Значения функций б/а, 92а) и #ска>. 


2а 6/23} 
1078 

т0-7 

1076 

1075 

10-4 

103 

0.0т 

0.т 

Т.О 5.844 
2.0 5.63Т 
3.0 5.450 
4.0. 5.298 
5.0 5.175 
10.0 4.843 
25.0 4.38 
100.0 4.11 
103 

104 

105 


12 


6/4 


2233888888 


5.450 


Е 


1.67.1093 
1.67.10-8 
1.659-10-7 
1.568 .10`6 
0.668-1075 
-0.775.1073 
-0.513-ТО-Е 
-0.129 
-0.т86-Т07Г 
0.164.т10-Г 
0.324-т0Т 
0.403 то! 
0.441 ТОТ: 
0.444. т07Г 
0.296«т07Г 
о.т00-т0-: 
о.ттт-то-® 


0.112-1073 


4, 


4.44109 
4.45.10-8 
4.498107 
4.980-т0-6 
9.764-1075 
5.454-10-3 
0.286 
0.781 
0.475 
0.399 
0.380 
0.385 
0.40т 
0.493 
0.705 
0.905 
0.9896 
0.9989 


0.1125.10`4 1.000 


+ 9 жа 
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Данные численных расчетов (см.таблицу ТГ) согласуются с ре- 
зультатами асимптотического анализа выражений (25) - (27) при 
20 -— Ои 24а - ©.. Так, ши а > е> имеем, 
опуская члены порядка () Г‹ жа)*] : 


м - „9 „1-21 
| ва) = 2 ^ жа р ая 2 Ф(ха)=1 Ра’ 


Отметим, что последнее выражение (для А (ха)) совпадает с 
результатом работы [2Т] , полученным для малых значений дзета- 
потенциала (ом. [21] , формула (3.29)). $ 

При 2%а -> О0с точностью до членов [Фу (ха) ] имеем: 


{ (жа) = 1, 3 ска) = 2, 4, (ха) = 5 жа. 


На рис. Г показаны графики : 9 1 в функции 
аргумента 64 га. 
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